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Abstract 

O 

, 

Establishing CP violation in B^{B^) — > decays requires a measurement of polariza- 
tion of the final lepton pair, or a precise determination of the 5° — *• and — > 
^ ■ rates. We first argue that if the amplitudes of these decays are dominated by the scalar 
j> '. and pseudoscalar Higgs penguin diagrams, as happens e.g. in supersymmetry with large 
I tan/?, the CP asymmetries depend practically on only one CP violating phase. This 
^ • phase can be large, of the order of the CKM phase, leading to large CP asymmetries in 
■ ^ ' the T~^T~ decay channel of -B^(-Bj]) mesons, potentially measurable in BELLE or BaBar 
experiments. Secondly, we show that the existing TAUOLA r-lepton decay library supple- 
mented by its universal interface can efficiently be used to search for B^{B^) t^t~ 
decays, and to investigate how the CP asymmetry is refiected in realistic experimental 
observables. 



1 Introduction 



40 years after its discovery in the Kaon system, CP violation has also been firmly estab- 
lished in -B-physics in a series of high statistics experiments which enabled sufficiently 
precise measurements of the relevant observables. While the reported small difference 
between the time dependent CP asymmetries measured in B ^ J/ ipKs and in B ^ (pKs 
decays [1], if confirmed, would confiict with the Standard Model (SM) prediction, there 
is as yet no convincing signal of a contribution to CP violation from physics beyond the 
SM. Such a signal may eventually be provided by yet more precise measurements and 
joint analysis of CP asymmetries measured in different channels. Equally important is to 
look for CP violation in channels in which no (or negligibly small) effects violating CP 
are predicted by the SM. Observation of a non-zero effect in such processes would be an 
unambiguous signal of new physics contribution to CP violation. 

From many possible channels^ we consider in this letter fiavour changing decays of 
the neutral B mesons into lepton pairs, B^g — > . This channel has attracted a lot 
of attention since it is very sensitive to new physics which affects the fe-quark Yukawa 
couplings [2-5]. Approximate and full one-loop calculations in the supersymmetric exten- 
sion of the SM with large tan /5 (the ratio of the vacuum expectation values of the two 
Higgs doublets) [4-6] showed that in such a scenario one can expect truly spectacular 
enhancement of the rates of the decays B^^^ . 

Moreover, new physics can also lead to observable CP violation in these decays, which 
is not predicted by the SM. For example, CP violation could manifest itself through 
non-equal leptonic decay rates of the tagged B^{t) and B^{t) states {B^{t) and B^{t) 
are the states which at t = are tagged as B^ and 5°, respectively). If polarization of 
final state leptons can be determined, additional information on CP violation could be 
provided by non-equal T{B^{t) Ipl) and T{B%t) Ip^) [or T{B%t) Ip^) and 
r(i?0(t) ^ Ipl)] decay rates [7,8]. 

Theoretically leptonic decays are particularly clean as the only nonperturbative quan- 
tities, on which their rates depend, are the 5° meson decay constants -Fb^^. Moreover, 
Fb^ ^ cancel out in suitably defined CP asymmetries. On the other hand it is not clear 
whether the CP violation in these channels can be accessible experimentally. First of all 
these decays have not yet been seen: the best upper limit on the B^ /^^/^~ branching 
fraction comes at present from the BaBar experiment [9]: 

5r(5^ ^ /iV~) < 8.3 X 10"^ , (1) 

which improves on previous limits of BELLE (1.6 x 10"^ [10]) and CDF (1.5 x 10"'' [11]). 
The 5° branching fraction is bounded by [12] 

Sr(5° ^ fi+fx-) < 2.7 X 10"^ (2) 

resulting from the combination of the CDF (5.8 x 10""^ [11]) and DO (4.1 x 10"^ [13]) 
(all limits are at 90% CL). The limits (jlj and ^ are still about 2-3 orders of magnitude 

^For example, also very small CP violation is predicted by the SM in the charm sector. 
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above the corresponding predictions of the SM: Br{B^ 1-^^ 1^ ) ~ 1-3 x 10^^° and 
Br{Bg — > ~ 3.6 x 10~^ [14-16]. The main uncertainties of these predictions come 

from the decay constants Fb^ and Fg^ which are known up to a precision of ~ 15% [17]. 
If the decays B^^ — > //^//~ occur at the rates as predicted by the SM, their detection 
will become possible only at the LHC. New physics (like supersymmetry) can increase 
significantly their rates to a level that they can be observed at BaBar, BELLE or Tevatron 
in near future. However, as we will show, in this case the ratio of time integrated leptonic 
decay rates of the B^{t) and B^{t) states to is unlikely to deviate appreciably 

from unity. Polarization measurement also seems very difficult in the case of the 
channel [18]. In the t+t~ channel the situation is quite different: large CP violating 
effects can be expected, and r polarization measurement is possible. So far this channel 
has not been much explored experimentally on account of difficulties with identifying r 
leptons. As a result, practically no experimental limits on Br{B^g — > t~^t~) exist despite 
the fact that the corresponding rates are expected to be a factor ml/mj^ ~ 300 larger 
than the ones for decays into final states. 

The purpose of this paper is twofold. Firstly, we point out that in realistic scenarios of 
large tan (3 MSSM, where the rates of — > decays are significantly increased to 
a level measurable at the running BaBar and BELLE experiments, the CP asymmetry in 
the B^{B^) — > T~^T~ channel can be quite large and potentially measurable. Secondly, we 
identify realistic experimental observables in which the CP asymmetry can be detected. 
In addition to the ratio of time integrated leptonic B^{t) and B^{t) already mentioned, 
we consider two complementary observables: the tt^ energies from r — > tti/ decays, and 
the acoplanarity angle between the decay planes of the p mesons which originate from 
r —>■ pp. The former is sensitive to the longitudinal, while the latter to the transverse 
polarizations of r's coming from B^ and B^ decays. We show how the existing TAUOLA 
package [19-21] and its universal interface [22-24] can be used to search for these 
decays and for the CP violation, demonstrating that the necessary tools for full Monte 
Carlo simulations are reliable and ready for use in the B^{B^) t^t' physics. 



2 General formulae 

The effective Lagrangian describing B^{B^) — decays can succinctly be written as 

Aff = B%iPi{hs,d + a.,d7')^z + BlMKd + cis,d-t'')^i . (3) 

To simplify the notation the subscripts d and s referring to non-strange and strange B^ 
mesons, unless explicitly written, will be omitted. Hermiticity (CPT invariance) implies 

6 = 6*, a = -a* . (4) 

The amplitudes of S° decays into two helicity eigenstates read 

AL^{llll\Bl,) = MB{a + hP) , 

AR={llln\Bl,) = MB{a-h^) , (5) 
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where = (1 — 4mf/Af^)^/^. Similar formulae with a and b replaced by a and 6, respec- 
tively, give the amplitudes Al and Ar for the corresponding decays. Evidently, if 
both coefficients a and b are simultaneously nonzero, r(i?° — > /^/^J^) 7^ r(i?° — > /]^/]^) but 
this is not yet a signal of CP violation. CP is violated, for example, if 



(6) 
(7) 



because the initial and final states on both sides transform into each other under CP [7]. 
As there are no strong phases involved, this can occur only through the mixing of the 
and mesons. In the standard formalism [25] the state B^^^J^) (-^phys(^)) which at 
i = is a pure B^ (5°) evolves in time according to 

l^iysW)=5+WI^°) + V(t)|5°). 

Neglecting the difference of the decay widths of the two B^ mass eigenstates one finds 

AM 



9-it) 



t, 



2 

AM 



(8) 



where AM = M^o^ - Mso <^ Mb = {Mbo^ + Mbo)/2. The ratio p/q (calculable in the 
SM or its extensions) is given by 



p _ 




M* 


Q 


V^12 ^ 


\Mi2 


12 = 







1 Im — — 

2 M12 



(9) 



with Hu = M12 + f Pi 
state, which initially was a 5°, decays at time t into left-handed leptons is therefore given 
by 



{ItlllB^^ysit)) = 9+it) Al + -9-it) Al , 



and the rates of B^ 

WpL\Bl^yM)\' - 



Ipl and 50 



IrIji decays are proportional to 



(10) 



1 + 



qg-it) 



P9+{t) 
P9-{t) 



At 



Q9+{t) 



Al 



Ar' 



A 



R 



\pg+{t)ALj 

+ 2Ref^^4^^ 
\qg+it) Ar^ 



The matrix elements for B^ 
by interchanging L ^ R. 



I'rIr and B^ Ij^lj^ can be obtained from the ones above 
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Since \ Al\ = I-^rI, I-^rI = \ Al\ and, as follows from ^ and 

Al \aJ ' 

one sees that the CP is violated if either 

^1 or Im(AL)^Im(A^^) , 



(11) 



where 



Al = — -r- , and Xr = — — . (12) 
pAl pAr 

The simplest quantities measuring the amount of CP violation are the asymmetries 
constructed out of time integrated polarized decay rates 

,1 _ dt r(i?pV3(t) ^ itii) - 1,'^ dt T{Bl^,,{t) ili-^) 

~ a: dt TiB^sit) - IPl) + III dt TiB^sit) - IP~r) ' ^ ^ 
,2 .._ ll:dtT{BXAt)^lpR)-ll:dtT{Bl^^M^lp-L) .... 

~ a: dt T{BXM - iPr) + III dt T{BXM - ipi) ■ ^ ^ 

and the ratio of integrated unpolarized decay rates 

The time interval (ti, ^2) can be chosen according to the experimental convenience. If the 
number of tagged events is not very large, or there is a large uncertainty in experimental 
determination of the decay time t, one can exploit the asymmetries A}.p = yl^p(0, 00) 
and Aqp = Aqp{0,oo) and the ratio Ri = -R/(0, 00) for which it is straightforward to 
obtain [7, 26] 



.1 2- l-R IJ - — -R J 





- \xp\ 


) - X Im (Al - X-p) 




2 + X2 + 


(|Ail + |A^^|)-xIm(AL + A^^^) 
-IAZ^l)-xIm(A^-AZ^) 


2 + x2 + ia;2 


(|A|| + |AZ^| 


) - a; Im (Ar + A^') 



and 



^ _ (lAP + l^i^n (1 + + ix2|gn - x{\AL\'lm{XL) + |^i,pIm(AH)} 

(|A|2+|^^|2)(l + ix2 + |x2|£|2)-x{|A|2Im(A^i) + |^«|2Im(Azi)} ' ^^^^ 
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where x = AM/T. If \q/p\ = 1, the relation (fTT| imphes \Xl\ = IXr'I- Moreover, A^ + A^^ 
is then real, and the formulae (fT6|) . (fT7|l and (fT8|) simplify to 

2 X Im Al 2 ~2 X ImAfl 



"^^P 2 + x2 + x2 |Az.|2 ' 2 + x2 + a;2 jA^P ' 

P _ (I + IArI') (1 + x^) - X {| APlm(A^) + l^^pim(A^)} 



The asymmetries A^P' ^cpi functions of Xl,r, are bounded from above by [7] 



Ah2 



< • (20) 



Since Xg > 20.6 for the B^-B^ system, the CP asymmetries in the leptonic Bg{B'^) decays 
can reach at best ~ 4.5%. In contrast, for the B^-B^ system, for which Xd = 0.771 ±0.012 
they can be as large as ~ 60%. This is fortunate, since B^B^ are copiously produced 
at BaBar and BELLE in a relatively clean environment (compared to the 5° production 
at hadron colliders). For this reason we will consider only the CP asymmetries in the 
50(5°) ^ /+/- decays. 

The quantities (fT3|l and (fT^ depend on asymmetries of B^ and B^ decays into longi- 
tudinally polarized leptons. In the case of the t~^t~ decay mode they are best identified 
by measuring the 71"^ energy spectra from r ^ vrz/ decays [27]. The density matrix for- 
malism outlined in Section |3| allows to construct also observables sensitive to transverse 
polarization of the final state r's [28]. These observables will prove complementary since 
in some scenarios the signal of CP violation can clearly be visible in the latter observables 
while hidden in the former (as Ri can be expressed in terms of polarized decay rates, Ri 
is then equal 1). 



3 Supersymmetry scenario 

If the Cabibbo-Kobayashi-Maskawa (CKM) matrix is the only source of both flavour and 
CP violation (as in the SM), then q/p ^ —V^lVtd{s)/ytbV*d{s) implying \q/p\ ^ 1. Moreover, 
Al,Ar oc V^lVtd{s), Al,Ar oc VtbVj^^g), so that Al and Xr are almost real. The CP 
asymmetries in the 5° decay are then negligible.^ In models of new physics one 

can have Im(A2,) ^ Im(A]j^) and/or \q/p\ ^ 1, but if the predicted rates of these decays 
are still of the same order as the SM predicts, the detection of the fi~^fi~ decay mode will 
become possible only at the LHC where most probably the measurement of the muon 
polarization will not be feasible. Detecting CP violation in this mode could be possible 
then only by measuring the ratio R^. The detection of the B^ — > r+r~ mode occurring at 
the SM model level {Br ~ 5 x 10^^) is possible at BELLE and BaBar ^ but the number of 

^In the SM they are estimated to be 0(10^^) [7] as a result of a small departure of \q/p\ from 1. 
^According to the SM prediction, among 5 x 10^ events of B'^B'^ pair production collected by these 
experiments so far there should already be some 50 events of decays into r+r^. 
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reconstructed events might be too small to detect any CP violation. On the other hand, 
at LHC where the number of — > t^t~ events will be larger, the identification of the 
r decay mode will probably be quite difficult and the CP asymmetry hard to detect. 

Much more promising situation can occur in the supersymmetric scenario with a large 
ratio of the vacuum expectation values of the two Higgs doublets, Vu/vd = tan/5 ~ 40-i-50. 
The coefficients a and b (and a and 6) in Q, and hence, the amplitudes of the B^^ ^ l^l^ 
decays can then receive important contributions from the Higgs penguin diagrams with 
s-channel and A° Higgs boson exchanges [4, 5, 8]. If the mass scale of the Higgs 
particles H'^ and (for Vu/vd 1 Mh ~ Ma) is not too high, of order < 500 GeV, 
the decay amplitudes can be dominated totally by these diagrams easily saturating the 
experimental limits (^,(121). This can happen even if the supersymmetric particle masses 
are quite large, say in the TeV range [5,29]. 

As illustrative examples of new physics we consider here two different supersymmetric 
scenarios: minimal (MFV) and non-minimal (NMFV) flavour violating, both with large 
ratio of VEVs. 

MFV Scenario: In this case no additional flavour violation in the sfermion mass matrices 
is assumed {i.e. the gluino-quark-squark vertices conserve flavour) but additional CP 
violation phases are introduced by complex parameters and At (the higgsino mass term 
and left-right top squark mixing, respectively). The relevant effective flavour violating 
couplings of the two heavy neutral Higgs bosons to the down-type quarks'' can be written 
as 



£f;-^"- = C VlVtd h 

+c v:^Vtb d 



Mw ' Mw fJ^* 



d 



b 



My/ fl Mw ^ 

where Pl,r = (1 =f 7^)/2 are chiral projectors and the coefficient C is given by 



9 .-„„2 n H2 



C = —^rTT- tan 3 -k 

4 M2 167r2 



with a dimensionless function of higgsino and stop masses if 2 of order C(l). The factor 
K oc tan (3 summarises some refinements in the calculation (resummation of tan /3 enhanced 
terms, for details see [6,29]). It is the factor tan^/? which makes these couplings so 
important. Combining these couplings with the and A^ couplings to and using 

{0\bP,d\B',) = -^Fs,-^ , (0|JP,6|5°) = +IFs,-^ , 

m^RABl) = +-Fb,^ , {0\dPnb\B',) = —Fb,^ , 

A Tflf) A Tflb 



^The effective flavour violating coupling of the lightest CP even Higgs to the down- type quarks is 
negligible. 
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one arrives at 



Ml / m, a; , m, AA _ ^, m. Ml A*, 



Ml \Mw Mwl^*J MvkMI ^ ' 

h-h*-n' T/*T/ / Al nid At\ _ , Ml A* 

Mfj \Mw /J' Mw /J' J Mw Mfj II 



with 



nil Fb ^ 3 ^ ^2 



16 M^Mwrub 167r2 



In this case the amplitude of the B^-B'^ mixing is not modified (in contrast to the one for 
B^-B^ mixing) [29], so that \q/p\ ~ 1 holds true as in the SM. One then finds 

AL~e A^~e (21) 

where the effective CP violating phase is given as 5cp = &xg{ii*Al). The time integrated 
asymmetries then read 

A^p = -0.09 X sin(25cp) , ^cp = -0-35 x sin(25cp) ■ (22) 

Since the sparticles giving rise to substantial Higgs penguin contributions to a and h (a 
and h) can be quite heavy, even order 0{1) phase of ^xAt needs not produce unacceptable 
electric dipole moments. 

NMFV Scenario: In this case squark mass matrices violate flavour and the corrections 
generating the flavour changing couplings of the neutral Higgs bosons ad A^ are 
dominated by gluino loops [5]. The relevant effective Lagrangian is then 



d 



with 



+D d fQ;d(^° + «^°)Pl + OiliH'' - a°)P J h , (23) 



m~ / rUb bd , M \ ^ rn- f nid , .db 



and 
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where is another dimensionless function of order C(l) of gluino and sbottom masses 
and K again denotes dominant higher order contributions [30]. The mass insertions S'^j^ 



7 



and etc. are the ratios of the diagonal entries of the down-type squark mass squared 
matrices to the average squark mass squared M|. One then finds 



a --a*- D' (a. + a*) ~ D'^^ ( "^S"" + '^S''"* 



b = b* = D' (a, - a*,) ^ D'i^'-i§ -^6^^ - 



with 



D = --gs9 ir^tan Pjj^-^ k — . 

The insertions 5^*^ and 5^ are not very tightly constrained. Typically the bound of 
order |5i'2.(i?_R) I < 0-2 TeV) arising from AM^^ is quoted [31]. This estimate, 

however, does not take into account the contribution of the so-called double-penguin 
diagrams [30,32-34], which can significantly affect the B^-B^ mixing amplitude. This 
contribution arises from the and A° exchanges between two effective (1-loop generated) 
vertices (f23| of which one vertex annihilates a right-chiral and the other a left-chiral 
quark^ and is therefore proportional to af,a^ oc S^^^d'^^ = ^il^/j/j (with the small d- 
quark mass neglected). Since this product is constrained much stronger [30], to avoid 
a potential confiict with the value of AM^^ and the time dependent CP asymmetry 
aj/^Ksi't) measured in i? ^ J/ipKs decay, we assume that either S''^j^ ^ 6^j^j^ or 6^^j^ <^ (J^"^. 
This leads to a ~ ±6 and, as in the previous scenario, to \q/p\ ~ 1. The asymmetries 
then read 



Al^p = -0.35 X sin(25cp) , A^p = -0.09 x sin(2(5cp) 



where now 



6cP = argiVtlVtd) - arg(m^/i5f^(^^)) . (24) 



It is interesting to note here that the CP asymmetries can be nonzero even if the phase 
of the CKM matrix remains the only source of CP violation {i.e. all supersymmetric 
parameters are real). Moreover, since |arg(V^^Vid)| is of order 1, the total phase violating 
CP needs not be small. 



In both scenarios, in which the ^ l^l^ amplitudes are dominated by the exchange 
of and Higgs bosons whose effective fiavour violating couplings to bd (or bs) differ 
only by a factor z, one gets 

a ^ b or a ^ ~b , 

^Contributions of H'^ and A° exchanges between vertices annihilating same chirahty quarks to the 
mixing amplitude are proportional to — « [32]. 



8 



(up to ^ 15%). For a = b the factors A^, and Xr defined in Eq. (fT^ simplify to: 

q a*l-f3 q a*l + (3 

Al = , Ar = ^ [2b 

p a 1 + p p a 1 — p 

and all CP-sensitive quantities depend on one effective phase which can be taken as 

Sep = -^arg(Ai) . (26) 

The immediate consequence of a = 6, with \q/p\ close to 1, is that for the /i^/i^ final 
state the parameters \Xl\ and |Ar| assume values ~ 4 x 10~^ and ~ 2.5 x 10^, respectively, 
since in this case /3 = (1 — 4m^/M^)^/^ is almost 1. As a result, the predicted asymmetries 
are very small: \Aqp\ < 2 x 10^^, I^cpI ~ 10"^. The same is true if a and b are somewhat 
split. In contrast, for the r+r^ final states, for which /3 = (1 — 4m^/M^)^/^ differs 
substantially from 1, we have \Xl\ ~ 0.15, \Xr\ ~ 6.7 and the maximal possible values of 
the asymmetries are 

(A^p)--| = 9% and |(v42p)"^^-| = 35% . (27) 

The comparison of magnitudes of possible CP violating effects in the ratio (fTSll for 
and T+r" decay modes is shown in Fig. [H where i?^ and Rr are plotted as functions 
of b/a for four different values of the phase Sqp (keeping arg(a) =arg(6) and \p/q\ = 1). 
The plots show, that the ratios Ri approach unity for a ~ ±b/3. (Vanishing of Ri for 
a ^ ±6/5 follows also from the formula (fT9|l if one takes into account that for a — > 6 ± /5 
Xr ~ l/(a =F b(3) whereas ~ |a =F &/?P-) Therefore, for a ^ ±b the deviation from 

unity of i?^ is tiny while for Rr it can be quite substantial. 

The asymmetries Aj^p (fT3|) . (fT^ and the ratio Ri (fT5|) are the observables in which 
the signal of CP violation (i.e. a non-zero phase 6cp) vanishes for a ~ ±bp. As the spin 
density matrix formalism will show, transverse polarization of the final state leptons is 
free from this problem. However, observables sensitive to the transverse polarization can 
experimentally be constructed only for r's which decay into hadrons. Thus, for a ^ ±b 
only the t~^t~ channel provides an interesting opportunity to look for CP asymmetry in 
the leptonic B'^{B^) decays. 

The coefficients a and b (a and b) are constrained by the experimental limit Eq. Jll), 
which in the case a ~ ±6 gives 

\a\ ^ |6| < 4.9 X 10^^ . (28) 

In the minimal flavour violation scenario, the parameters a and b are also constrained 
indirectly^ by the limit Q. Taking {mT-/m^){\Vtd\/\Vts\) ~ 4 this constraint is satisfled if 

\a\ ^ |6| < 2 X 10"^ . (29) 



^The limit imposed by the B^g-B^^ mass difference [29,32] can be avoided on account of a different tan (3 
and Ma^h dependence as compared to the the Higgs penguin diagram contributions to B'^{B^) l^l^ 
decays [16]. 



9 




Figure 1: The ratios and Rr as functions ofb/a for the phase 6cp = — |arg(Ai) = 0.1 
(solid line), 0.3 (dashed), 0.5 (dotted) and 0.75 (dash- dotted). Ri{—Scp) = Rl^{5cp). 

In our simulations in Section El we conservatively set \a\ = \b\ < 10~^ and treat both 
scenarios simultaneously, as all what matters are the values of \a\ = \b\ and the single 
CP violating phase 6cp. In the MFV case 6cp = arg(yU*A^) while in the NMFV case it is 
given by Eq. ^M : in both cases the phase can be of order 1. 

Finally let us notice that the SM as well as supersymmetric box and penguin 
contributions spoil the exact equality a = ±b. In addition, a finite difference of and 
masses also splits these two coefficients. It is therefore likely that \a\ and |6| and the phases 
of a and b differ from each other by some lO-i-15% even for \a\ and |6| saturating the bounds 
JH), Nevertheless, our simplified analysis will demonstrate that in supersymmetry 
there are good reasons to expect substantial CP violation in B^{B^) t^t~ decays. 
Therefore in the following sections we present tools which can be used to search for these 
decays and look for CP violation. 



4 Spin density matrix formalism 

To see how the CP asymmetry in — > r+r^ decays are refiected in realistic observ- 
ables we use the TAUOLA r-lepton decay library which allows us to simulate fully the effects 
of T polarization. The exhaustive description of the method and numerical algorithm is 
given in papers [19] and [22]. The input to the TAUOLA Universal Interface [24] is the 
spin density matrix of the t'^t~ system resulting from the decay of a neutral particle. In 
this section we collect the necessary formulae for this matrix for a r+r" pair originating 
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from B^B^). 

The time dependent S^'-meson mixing can easily be dealt with by introducing time 
dependent effective factors aes, bes, flefr and bes defined as: 

aefr(t) = a g+{t) + a - g^{t) = a g+{t) - a* - g^{t) , 
p p 

b,s{t) = b g+{t) + b^g4t)=b g+{t) + b* ^ g_{t) , 
p p 

_ _ p p 

a^sit) = a g+(t) + a - g-{t) = -a* g+{t) + a - g^{t) , 
q q 

Ks{t) = b g+{t) +b^g4t)= b* g+{t) +b^g4t), 
q q 

so that the instantaneous -B° widths into left- or right-handed r's read 



r(5ShysW - rlr,) = ^(3 Mt) + (3 b^^' , (30) 



IGvr 



nS^M - r^r^) = laefr(t) - /3 b,^{t)\' , (31) 

and those for S° are given by similar formulae with a^sit), b^sit) replaced by aes{t), bcsit)- 
CP is violated because in general aes(t) ^ — a*g(t), and 6efr(i) 7^ ^cff(^)- 

The spin weight for the r+r^ pair originating from B^ decay at time t is given by 

-t = i^^^^i^, (32) 

where 

T(si,Pi,S2,P2) = Tr P(si) (^1 +m/) (6cfT + aeff7^)P(s2) {i)2 - m) (&eff -OcffT^) 

with -P(s) = |(1 + 7^/), and pi, Si (p2, -52) are the momentum and spin four-vectors of 
the T~ (r+) lepton, respectively. The r-lepton spin four-vectors Sa in the B^ (B^) rest 
frame are related to the spin three-vectors (Tq in the Ta-lepton rest frames as follows^ 

3? = <yl , si = -^(3 , 

^ 2mi ^ ' ^ 2mi ^ ' 

Mr 

Combining (p^ with the r-lepton decay matrix elements, the r-lepton rest frame spin 
vectors ai and (J2 get replaced by the polarimetric vectors hi and /12 which are determined 
solely by the dynamic of the considered r decay process. 

^In the rest frame of the z-axis is ahgned with r~ momentum, exactly as in TAUOLA Universal 

Interface. However, here as a particle number 1 the r~ lepton is taken (not t+), resulting in the 
transposition of the matrix R^^ in the TAUOLA Universal Interface. 
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The spin weight for the complete event [B ^ tt ^ decay products) can be written 
in the form 

^T=-il+ E Rr,Khi+ E RroK+ E R^^hi] . (34) 

\ i=x,y,z j=x,y,z i=x,y,z j=3:,y,z 

Expanding (p!^ and comparing with we find 

RqO = +1, RxO = RyO = Rox = Roy = 0, 

Rzz 1; Rxz Ryz Rzx Rzy 0, 

_ _ 2Re(aeff&:ff) P 

-rtOz — —-n-zO — 77 12,32 1 I 12' V"^"^/ 

P __D 2Im(aeff&:ff) D _p _ \beS?P^ - 

rLxy — rLyx — I 12/32 I U 12' ^xx — -ttyy — \2iQ2j_\n |2 ' 

Eq. (|3i|l with i?^,^ replaced by computed as above but with aesit) and 6eff(^) replacing 
aes{t) and &efr(^), respectively, gives the spin weight for the events from -B° decays. CP 
violating effects in the t^t~ and t^t~ decays are absent if Rqz = Rzo, 

Rxy Ryx and Rxx Rxx ■ 

For simulations of the time integrated measurements, the time averaged matrix {R^u) 
has to be used ^ 



{Rfiu) 



-■phys ^ 

and {Rfj,u) given by a similar formula. The asymmetry (fT3|) is then given by 

,1 _ (1 - {Rz0))Tint - (1 + {Rz0))Tint . . 

(i-(i2,o))rint + (i + (i?.o))rint ' ^ ' 

where 



int 



J dt T{Bl^^,{t) r+r-), T,^, = J dt r(5°,y,(t) ^ r+r^ 



^cp(^i'^2) defined in (fT^ is given by (f36|) reversing the signs in the brackets. It is 
also easy to check that for \q/p\ = 1 and a = ±b one has {Rzo)Tint = — (-Rzo)fmt 
Re(aeff&efr) = —^^{(iesKff)) ■ Then the two factors: (-R2o)rint and {Rzo)Tmt cancel out in 
the numerator of (f36|l and of the similar formula for AQp(ti,t2). As a result, nonzero 
asymmetries A^p are possible only if Tint 7^ Tint, which in view of the equalities |aefr(t)| = 
|6cfr(^)|, |^cfr(^)| = |aefr(t)| requires P significantly different from 1. This again confirms our 
observation made in Section [21 that for a ~ ±6 the asymmetry in the fi~^fi~ channel is 
suppressed. 

* Simulations for time non-integrated measurements, with time-dependent decays explicitly generated 
are also possible. 
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Since in the limit a = ±bf3 the two integrated rates Tint and Tint are equal one gets 
{Rzo) = —{Rzo), which means that even for 5cp 7^ there can be no CP violating effects in 
the observables sensitive to the longitudinal polarization of r's (nor in Rr). In contrast, it 
is easy to check by using the explicit analytical expressions, that in this limit the elements 
XX and xy of these matrices need not satisfy (Rxy) = ~{Rxy) and (Rxx) = {Rxx)- Hence, 
the observables sensitive to transverse r polarization can reveal CP violation even if the 
observables introduced in Section [2l fail to signal it. 



5 Results of the Monte Carlo simulations 

As 5° T^T^ decays have not yet been seen it might seem premature to study differ- 
ential distributions in this channel, including r polarization. Nevertheless, our analysis 
can serve as a good starting point for the future experimental work if indeed accumulated 
experimental samples would turn to be large enough. By providing a Monte Carlo tool 
useful in calculating e.g. detector responses, our study may also be considered as bench- 
marks for the simulations to be used in setting the upper limit for —>■ t^t^ branching 
ratios. 

There are some similarities between this study and the one aiming at assessing mea- 
surability of the Higgs boson parity at future Linear Colliders [28]. The difference is that 
now the coefficients a and b (and also a and b) in Eq. can be complex, while in [28] 
they were taken to be purely imaginary and real, respectively. As a result, in the formula 
Eq. terms linear in the polarimetric vector components also appear. Therefore we 
have extended the TAUOLA Universal Interface to include such a possibility as well. 
The algorithm for simulating decay into r leptons is nearly identical to the one for 
decay. Changes necessary for implementation in the Universal Interface of the TAUOLA 
Monte Carlo library are limited to the replacement of Higgs identifier with the one for B^ 
and B^ mesons, and the values of the spin density matrix with the ones computed in the 
preceding section. 

The formulae given in the previous section are general. However, motivated by the 
two supersymmetric scenarios discussed in Section O in our numerical study we will show 
first the results obtained in the limit a = ±b. Then we will discuss effects of a small 
departure from this relation. 

We will limit ourselves to two observables, which are known to provide valuable and 
complementary information on the spin state of decaying r lepton pairs. 

TT^ energies: As the first observable we take the 7r+ and n~ energy spectra in the decay 
channels r+ tt+z/^ (or n^Ur). Since they refiect the longitudinal polarization of 

the individual leptons, the spectra are sensitive to Rzo and Rqz as can be inferred from 
the expression (jSH), i-e. they are sensitive to Re(aefr6*fj) as follows from (jHSI)- The CP 
violation is refiected in the difference between the energy spectrum of 7r~(7r+) originating 
from B^{B^) and the energy spectrum of 7r"*'(7r") originating from B^{B^). This observable 
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was exploited before, for example in the /'^* t^t~ study (for a review of the method 
and its extensions, see e.g. Ref. [35]). We will show energy spectra in the rest frame of the 
B^{B^) meson assuming that the reconstruction of the event kinematics in the BELLE 
and BaBar experiments is sufficiently good for that purpose, that is, that the momenta 
of r decay products in the rest frame of _B°(_B°) can be reconstructed with the precision 
of a fraction of the r mass. 

Acoplanarity angle v^*: As the second observable we use the acoplanarity angle (p* 
between two planes spanned by the momenta of decay products of —>■ tt'^tt^ coming 
from decays of both r leptons into pUr [28]. This quantity is sensitive to correlations 
between transverse components of r-lepton spins (i.e. to the elements Rxx and Rxy which 
in turn probe Im(aeff6*fj), as can be seen from (|35|l ). For the definition of the acoplanarity 
the orientation of decay planes and pion momenta has to be properly taken into account. 
The acoplanarity angle ip* is defined with the help of two vectors n± normal to the planes 
determined by the momenta of pions which originate from p^ decays: ri-t = p^ri x p^o. 
If cosf = I """"I'l"" then 

^ |n+||n_ 

* ^ f ^ for sgn(p^- ■ n+) < 

^ [ 27r - ^ for sgn(p^- ■ n+) > ^ ' 

making the full range of the variable 27r of physical interest. Note that under CP, 

(f* ^ 271 — if* since the condition sgn(p^- ■ n+) is always evaluated from the orientation 
of 7r~ momentum with respect to the vector n+. In addition, we also have to sort events 
depending whether yiy2 > or yiy2 < 0, where 

Et^+ + E^^o E^- + Ej^o 

since otherwise the spin correlations are washed out, as explained in Ref. [36]. The best 
would be to use in (j38|l the energies of vr^ and vr^'s in the rest frames of the corresponding 
leptons, but they are not directly measurable. Since in the B^[B^) rest frame the 
r leptons are only mildly relativistic, the difference of pion energies in this frame and 
respective rest frames of should not be very important. The acoplanarity distribution 
is then evaluated in the rest frame of the p^ p~ pair, but with the energies of and 
7r°'s in (|38|) taken in the rest frame of the B^(B^). The CP violation is reflected in the 
difference between the distributions of the acoplanarity angle ip* measured in B^ decays 
and the angle 2n — ip* measured in B^ decays for the same signs of yiy2- 

Fig. [21 shows the pion energy spectra and the acoplanarity distributions assuming 
\q/p\ = I, a = b = 10~^ and the CP violating phase 6cf = 0.7. For all plots the same 
number of 5 x 10^ r+r~ events from 5° and B^ decays has been generated with TAUOLA, 
although for the parameters chosen the ratio Rr = 1.32, see Fig. [H In the upper left 
panel the energy spectra of vr" from B^ decays (thick line with the slope proportional to 
{Rzo) [37]) and of 7r+ from B^ (thin line; slope oc {Roz)) are shown, while in the lower 
left panel shown are the spectra of from 5° decays (thick line; slope oc {Rqz)) and of 
7r~ from B^ (thin line; slope oc {Rzo)). The harder vr" energy spectrum from 5° decays 
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Figure 2: Results for the CP violating phase 5cp — 0.7 and a — b. Left panels: Single tt^ 
energy spectra in B^{B^) — > r^T~ , tt^ In the upper (lower) panel the thick 

line corresponds to the energy spectrum of (ofn^) from decays and the thin line to 
the energy spectrum ofii^ (of7T~) from B^ . Spectra are plotted in the rest frame of B^{B^). 
Right panels: acoplanarity distributions of the p~ decay products in B^[B^) t'^t~ , 
P^Vt{Pt), — > 7r^7r°. The thick lines correspond to the acoplanarity angle </?* 
measured in B^ decays and the thin ones are for the angle 2i\ — 99* measured in B^ decays. 
The acoplanarity angles are defined in the rest frame of the p~ pair. Events in the 
upper (lower) panel have yiy2 > (yiy2 < 0). 

than 7r~^ from 5° (i.e. larger slope of the thick line) in the upper left panel indicates 
that Br{B^ — t^t^) > Br{B^ '^l^l)^ which is a clear signal of CP violation. In the 
acoplanarity plots (right panels) thick lines correspond to the distribution of (p* measured 
in B^ decays, and the thin lines to the distribution of 2n — (p* measured in S° decays; in 
the upper right panel yiy2 > 0, and yiy2 < in the lower right one. The shapes of the 
thick and thin lines are described by the formulae 

Nb{p>*) = const - sgn(yiy2)^iiC0s((/?* - Sr) for S° t^t" 
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where Ar = {R^ 



const — sgn(yiy2)^R cos(27r — if* + 5 



R) 



for S° 



r r 



_|_ 132 



)^/^, sin(5i? = Ra^y/AR, cos5r = Rxx/Ar and Aj^ and 6r are given 
by analogous formulae with Rij replaced by the Rij. Different shapes of thick and thin 
lines seen in the right panels of Fig. [2l again indicate CP violation. In both energy and 
acoplanarity plots the CP violation is clearly seen and should be measurable even for 
small statistics. 

Note also that if upper and lower plots are combined {i.e. no sorting according to the 
pion charge or sign of yiy2 is made), all CP asymmetries are lost. Since the lower plots 
are simple reflections of the upper ones, in the following only plots corresponding to the 
upper panels of Fig. [21 are shown. 
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Figure 3: ^45 in the upper panels of figure\^ but for Sqp = 0.3 and a = b. 



Fig. ini shows the corresponding distributions for the CP violating phase Scf = 0.3 
keeping still a = b and \q/p\ = 1. It is apparent that with decreasing \6cp\ the signal of CP 
violation deteriorates (especially in the pion spectra) and the possibility of distinguishing 
pion spectra and acoplanarity distributions from and 5°, and hence the CP violation, 
would require increasingly large statistics which may not be attainable at BELLE and 
BaBar without major upgrades. 

As we discussed, the relation a = ±6 is only approximate. Therefore, in Figure 21 (0) 
we show the n"^ spectra and acoplanarities for a greater (smaller) than 6, but keeping as 
previously their phases equal (and assuming \q/p\ = 1). We take 6 = 0.8 a in Fig. [Hand 
a = 0.8 b in Fig. O In both flgures the single CP violating phase 6cp = 0.7. 

Fig m shows that for b = 0.8 a with the same value of Sep the CP violating effects in 
TT^ energy spectra are enhanced, while in the acoplanarities they are only slightly affected 
compared to the case a = b, c.f. the upper right panel of Fig. (21 

On the other hand, for a approaching bp (for B t~^t~ decays /? ~ 0.74) the effects 
of CP violation in the tt^ energy spectra disappear. This is clearly seen by comparing the 
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right panel of Fig. El with the upper left one of Fig. 01 This agrees with our observations 
following Eqs. ([25|l and (f26|l and with the discussion in Section (HI In contrast, the acopla- 
narities shown in the right panel of Fig. (HI clearly indicate the CP violation even for a ~ 6/3 
confirming our discussion in Section IH This clearly demonstrates the complementarity of 
the energy and acoplanarity distributions as a means to detect CP violation. 
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Figure 4: As in the upper panels of figure\E but for 6cp = 0.7 and b = 0.8 a. 
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Figure 5: As in the upper panels of figure\B but for 6cp = 0.7 and a = 0.8 b. 



6 Conclusions 



In this letter we have investigated possible signals of CP violation in the decays _B°(_B°) — > 
r+r~. We have developed the necessary formalism and numerical tools allowing to apply 
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the TAUOLA r-lepton decay library together with its universal interface to simulate 
fully the effects of the polarization of r"*" and r~ originating from such decays. 

We have argued that in the interesting new physics scenario of supersymmetry with 
tan/3 ~ 40 50, in which the rates of B^{B^) —>■ t^t~ decays are enhanced and could be 
detectable in the SLAC and KEK -B-factories, the dependence of the CP asymmetries on 
the model parameters simplifies. Moreover, the CP violating phase needs not be small. 
In the non-minimal flavour violation case it can be of the same order as the phase of the 
Vtd element of the CKM matrix. Therefore the CP asymmetries can be quite large as 
opposed to the B^{B^) — > decays in which they are kinematically suppressed. 

By using Monte-Carlo simulations we have investigated the possible effects of CP 
violation in two realistic experimental observables and demonstrated that they might be 
detectable if the CP violating phase is reasonably large, i.e. 0{1). 

Since the decays B^{B^) — > r+r~ have not been discovered yet, we have not discussed 
the statistics requirements nor attempted at including in our analysis any systematic or 
detector uncertainties. It is clear that once these decays are discovered, other r decay 
channels than the ones investigated here can be analysed jointly to give additional in- 
formation on the polarization of r's. Our numerical tools are prepared for that. The 
tools can also be applied to determine the upper limits on the branching fraction of the 
B^{B^) — > r+r" decays by the BaBar and BELLE collaborations. 

As a final remark, we point that our analysis can be taken over to Higgs boson pro- 
duction at linear colliders with its subsequent decay to r pairs, which as yet has not been 
simulated in connection with the complex scalar and pseudoscalar couplings. 
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